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Abstract 



A quantum kinetic equation has been derived for the description of pair pro- 
k> , duction in a time-dependent homogeneous electric field E{t). As a source term, 

H I the Schwinger mechanism for particle creation is incorporated. Possible particle 

production due to collisions and collisional damping are neglected. The main 
result is a closed kinetic equation of the non-Markovian type. In the low density 
approximation, the source term is reduced to the leading part of the well known 
Schwinger formula for the probability of pair creation. We compare the formula 
obtained with other approaches and discuss the differences. 



1 Introduction 

Particle production in ultrarelativistic heavy-ion collisions raises a number of challenging 
problems. One of those interesting questions is how to incorporate the mechanism of 
particle creation in a dynamical kinetic theory [|^, |^, |^. In the framework of a flux 
tube model, a lot of promising researches has been carried out. In the scenario where a 
chromo-electric field is generated by a nucleus-nucleus collision, the production of pairs 
can be described by the Schwinger mechanism [^, ^ |[. The charged particles produced 
generate a field which, in its turn, influences the initial electric field and may cause plasma 
oscillations. Recently the interesting question of this back reaction has been analyzed 
within a field theoretical approach 0, ||, |]. The results of a simple phenomenological 
consideration based on kinetic equations and the full field-theoretical treatment |T0, llll 



12, 13| agree with each other. The source term which occurs in such a modified Boltzmann 



equation was derived phenomenologically in [|I^ . However, the appearence of this source 



term in relativistic transport theory which may be of non-Markovian character [|1^, |T^ is 
not yet fully understood. 

In the present work, a kinetic equation is derived in a dynamically consistent way 
for the time evolution of the pair creation in a time-dependent and spatially homo- 
geneous electric field, a process referred often to as the Schwinger mechanism. This 
derivation is based on the Bogoliubov transformation for field operators between in-state 
and instantaneous states. In contrast with phenomenological approaches, the source 
term of the particle production is of non-Markovian character. The kinetic equation 
derived reproduces the Schwinger result in the asymptotic limit and in the low density 
approximation. 

The paper is organized as follows. In Section 2 we sketch some important features of 
our approach to be used in Section 3 for deriving the basic kinetic equation describing 
particle creation in strong electric fields. Finally we discuss the properties of the equation 
obtained. 



2 Model description 

We use a simple field-theoretical model to treat fermions in an external electric field. 
To be described by the vector potential A'^ = {0,0,0, A^(t)) [] the electric field E^(t) is 
assumed to be time-dependent but homogeneous in space and hence E^ = E'^ = 0, 

E{t) = E%t) = -A^{t) = -dA^{t)/dt . (1) 

This quasi-classical electric field interacts with a spinor field ip of fermions. We look 
for solutions of the Dirac equation where eigenstates with the positive and negative 
frequencies are represented in the form f\ |T^, |T8|, ^ ^, ^ : 



4fi^) = L-'/' 



z7°9o + 7Vfc-e7^A3(t) + m 



X^^\p,t) Rr e^^^ 



^We use the units h — c= I and the metric is chosen to be g'^''' — diag(l, —1, —1, —1). 
^It is convenient to work in the discrete momentum representation assuming the limited volume of 
the system, V = L^. 



where /c = 1, 2, 3 and Rr (r = 1, 2) is an eigenvector of the matrix 7*^7 
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so that R^Rs = 26rs ■ The functions x (p^t) are related to the oscillator-type equation 



x^^Kp,t) 



uj'{p,t)+teA3{t)) x^^KP,t) 



(4) 



with uj'^ip^t) = /i^ + Pg^, /i^ = m? + p\ being the transverse mass and -P3 = Ps — eA^{t). 
The solutions x^^\P) ^) of Eq.(^) for positive and negative frequencies are defined by their 
asymptotic behavior at t — > — 00 






exp ( ± iu^ (p) t) 



(5) 



where oj^{p) = lim uj{p,t). In the following we will sometimes drop the momentum 

t— » — CXD 

dependence of the frequency and use the short notation uj{p,t) = u}{t). Note that the 
system of the spinor functions @ is complete and orthonormalized. So, the field operators 
ip{x) and ■?/'(x) can be decomposed in the spinor functions (0) as follows: 



'?/'(x) 
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The operators ap^.' and apr describe the creation and annihilation of electrons and 
positrons in the in-state at t— *> — 00 with the vacuum state \Oin>- The evolution will affect 
the vacuum state and mix the positive and negative energies resulting in non-diagonal 
terms of eigenstates responsible for pair creation. The diagonalization of the Hamiltonian 
for a Dirac-particle in the homogeneous electric field (0) is achieved by means of the 
time-dependent Bogoliubov transformation 
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with the condition 



1 . 



\ap{t)\' + \f3p{t)\' 

Here, the operators bpr and bpr describe the creation and annihilation of quasiparticles at 
the time t with the instantaneous vacuum |0f >. The substitution of Eqs.(0) into Eqs.® 
leads to the new representation of the field operators, e.g. 



ipi^x) 



E 

r,p 



,i-), 



^ pr \X) Opr [Z ) ~r ^ pr \X ) 0_pr\Jyj 



pr 
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(9) 



'Below we follow the notation of Ref. |2S 



The new basis functions '^pr (x) are connected with the former ones (0) by a canonical 
transformation hke Eqs.(^ 



'ippr\x) 



■,(+) 



r(-) 



«p(t)^^;^-/3|(t)^ 
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^P^Vix) = a*Jt) ^'p^> + Ppit) ^ 



pr 



(10) 



Therefore one may assume that the functions \I/^r have a spin structure similar to tjj, 
given by Eq.(D, 



■'pr 



¥^\x) = L-'/' 



n°5o + 7Vfc - e7%(t) + m 



(,(±), 



X 



Rr e^'^^'^e'P^, 
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where the dynamical phase is defined as 

e{t) = f dt'ujit') 

Jtn 



(12) 



and (j)p are yet unknown functions. Substitution of ([TT|) into Eqs.(|TOD leads to the 
relations 



X^-\p,t) = a,(t) 0H(t) e-ew - /3|(t)0(+)(t)e^eW, 
xW(p,t) = a;-(t)4+^(t)e^®W + /5p-(t) 4-)(t) e-«W 



(13) 



Taking into account that the functions x (Pi't) are defined by Eq.(^), we can find the 
equations for functions ap{t) and Pp{t). According to the Lagrange method it is possible 
to introduce additional conditions to Eqs.([T3|) 



x^-\p,t) 



X^^\p.t) 



-iuj{t) 



iuj{t) 






(14) 



Let us differentiate these equations repeatedly and make use of Eqs.(§) and ([131) . Then, 
choosing the functions 4 (t) in the form 



^?\t) 



Lu{t)±Ps 



we obtain the following equations |18|, |19|, ^ pl[] : 

ap{t) = 

m) = 



(15) 



^^('^)^^ o*(.^ 2ie{t) 
eE{t)n 



(16) 
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One should note also the usefuU relation for these transformation coefficients 
2^u;it) (uit) - P3) ^ ^ 

It is convenient to introduce new operators 

cg^(t) = 6g)(t) exp(±ze(t)), 



(17) 



It is easy to show (e.g., see p2[) that these operators are consistent with the rules of the 
Hermitian conjugation for creation and annihilation operators, i.e. 

[cg^(t)]+ = cWW (19) 

and satisfy the Heisenberg-like equations of motion 






(20) 



where H(t) is the Hamiltonian of the quasiparticle system 

H{t) = '£u{t)(c^\t) cH(t) - cL-,l(t) 4i(t) ) . (21) 



r,p 



The first term on the r.h.s of Eqs.(pOD is caused by the unitary non-equivalence of the 
in-representation and the quasiparticle one. The commutation relations 

are fulfilled for the operators ([TsD . 

3 Dynamics of pair creation for fermions 

Now we shall consider the kinetic equation for the distribution function of electrons with 
the momentum p and spin r 

Mp,t) =< Oin\bi^\t) 4;^(t)|0,„ > = < Oin\ci'r\t) 4;^(t)|0,„ > . (23) 

According to the charge conservation the distribution functions for electrons and positrons 
are equal fr{,P,t) = fr{p,t), where 

fr{P,t) =< Ojb'^^lit) 6L"pU^)|0m > = < 0,n\c^-l{t) C^Sl{t)\Oin > • (24) 



The distribution functions ( PB| ) and (|2^) are normalized to the total particle number 
N{t) of the system at the given time moment, 



Y.friP,t)=J2friP.t)=Nit). 



(25) 



r,p 



r,p 



By the time differentiation of Eq. (p3D we get 
dfr{p,t) enE{t) 



dt 



2uj\t) 



^^+\p,t) + ^[-\p,t) 



(26) 



We have used here the equation of motion (^) and evaluated the occurring commutator. 
The functions ^)^' in Eq.(|26|) describe the creation (^^ ) ^^"^ annihilation ($}."■') of an 
electron-positron pair in the external electric field E{t). They are given as 



<l>(+)(p,t) = <0.n\*4t\t)c^_^l{t)\0,n> 



(27) 



It is straightforward to evaluate the derivatives of these functions. Then, by applying the 
equations of motion (pO]) , we get 



c?$(^)(p,t) _ efiEjt) 
Jt ~ 2cj2(t) 



2fr{p,t)-l 



±2iu;{t) <^i'^\p,t) 



(2J 



where due to charge neutrality of the system, the relation fr{p,t) = fr{,P,t) is used. The 
solution of Eq. (^) may be written in the following integral form: 



*^^np,t)=f r dt 



>eE{f) 



uj\t 



2fr{P,t') 



,±2i[e(t)-0(t')] 



(29) 



Here we proceeded to the limit to ~* "Oo for functions &{t) and B(t') in r.h.s. of Eq.(p9|) 
(see the definition (|T^)) assuming that lim A^(t) = A^ = 0, i.e. the fields $^^'*(p, t) 

vanish at t — ^ — cxd. Inserting Eq.(p9|) into the r.h.s of Eq.(^) we obtain 



dfr{p,t) e^iE{t) ft ^^,efiE{t') 



dt 



2uj\t) 



LU^{t 



l~2friP,t'] 



COS 2[0(t)-e(t')] . (30) 



It is seen from (|30|), that the distribution function does not depend on spin, so fr = f ■ 

Now let the volume of the system goes to infinity, V^ ^ oo. To do this limiting 
transition, we use the rule 

(|)E...-/^v.., 

take into account the normalization conditions (|25|) and introduce the distribution func- 
tion -F(p, t) normalized to the particles number density, i.e. 



{2n)-'g d'pF{p,t)=n{t). 



(31) 



where g = 2 is the degeneracy factor for electrons. Then, with the substitution f{p,t) -^ 
F{p, t), and after replacing the canonical momentum p by the kinetic one, p -^ P = p~eA, 
where the 3-momentum is now defined as P{p^^p'^^ P^), the kinetic equation (^) is reduced 
to the final form: 



dF{P,t) dF{P,t) 



dt 
with the Schwinger source term 



dt 



+ eE{t) 



dF{P,t) 
dP. 



SiP,t) 



S{P,t) 



2u^{t) 



dt 



,e^iE{t') 



l-2F{P,t'] 



cos 2[e(t)-e(t')] 



(32) 



(33) 



Recently, a kinetic equation similar to (0) has been derived within a projection operator 
formalism for the case of a time-independent electric field in Ref. [|T^ where it was first 
noted that this source term has non-Markovian character. The presence of the Pauli 
blocking factor [1 — 2F{P,t)] in the source term has been obtained earlier in Ref. p. 
We would like to emphasize the closed form of the kinetic equation in the present work 



where the source term does not include the anomalous distribution functions (P7|) for 
fermion-antifermion pair creation (annihilation). 



4 Concluding remarks 

Based on microscopic dynamics, the closed kinetic equation ( P^D has been derived for 
fermion evolution in a time-dependent, homogeneous electric field with the Schwinger 
source term (^) to be characterized by the following features: The kinetic equation (^) is 
of non-Markovian type due to the explicit dependence of the source term on the whole pre- 
history accounting for a memory effect. The difference of dynamical phases, 0(t) — Q{t'), 
under the integral (|30| ) generates high frequency oscillations. The appearence of such a 
source term violates the time reversal symmetry and may lead to entropy production due 
to the pair creation. As was noted in Refs. [^, |T^ , this may result also in oscillations of 
the relevant entropy and thus in temporaryn (on the scale of the memory time) violations 
of the H-theorem. The phase space occupation of the created fermions is taken into 
account by the Pauli blocking factor 1 — 2F{P,t') in (|33|). 

It is worthwhile to compare our result with the equation obtained phenomenologically 
in Refs. [|, 



S^^{P,t) = -[1 - 2F(P,t)] \eE{t)\ In 



1 — exp 



\eEit)\J 



S{P3 



(34) 



By comparing Eqs. (|33|) and (Q) we see that the essential difference between them is 
in the non-Markovian character of Eq. (|33[) . This rather intriguing feature leads to an 
accumulation of created particles on the background of rapid oscillations. 

Finally let us check the semi-classical limit for the case of a constant electric field, 
E{t) = E. From Eq . (|33D at t ^ -|-oo and in the low density limit F{P,t) <^ 1, the pair 
production rate is reduced to the following formula: 



S^^ = lim (27r)-^g f d^P S(P,t) = ^^exp 



JeEl 



(35) 



which covers the well-known Schwinger result obtained in the semi-classical approxima- 
tion. It is of interest to note that in contrast with this semi-classical limit the source term 
(0) is not positively defined, in general. 

Our consideration can also be applied to scalar systems described by the Klein-Gordon 
equation. Under the same assumptions we obtain Eq. ( ^2]) (see Appendix A ) with the 
source term 

^(^■') = flr/l*'^ [i + 2F(p,0]™(2|ew-e(0l) . (36) 

The difference between the sources ( ^31) and (^) is only in the phase space occupation 
factors due to the different quantum statistics of particles involved: the Pauli blocking 
factor 1 — 2F{P,t) which is present in the source term for the creation of a fermion pair 
has to be replaced by a Bose enhancement factor 1 + 2F{P, t) when a pair of bosons is 
created (stimulated pair production). 

Besides of the numerical solution of the derived kinetic equation, in our future studies 
we would like to incorporate collisions in this formalism and apply our method for the back 
reaction problem. In this case the time-dependent electric field must be defined by means 
of the Maxwell's equation with the spatially homogeneous current j{t) = —{l/A7T)E(t). 
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Appendix A 

Here some details are given to derive of the kinetic equation with the source ( PE| ) describing 
the vacuum creation of scalar mesons in a strong electric field. 

The solution of the Klein-Gordon equation in the presence of the electric field defined 
by the vector potential A^ = (0, 0, 0, A^{t)) is taken in the form [O, WB, M, 120, ^^ 



Wr^^ = r-3/2r9,,, (r,\]-y2 P^^Pn(±)rr, 



[X 



where the functions g'^^\p^ t) satisfy the oscillator-type equation with a variable frequency 

g'^^\p,t) +uj\p,t) g'^^Xp^t) =^ . {A2) 

Solutions of Eq.(A2) for positive and negative frequencies are defined by their asymptotic 
behaviour at t — > — oo similarly to Eq.(^. The mesonic functions (Al) are normalized 
(in the discrete momentum representation) in the following way: 



I 

V 



J c?3x0f )*(x) do (j)f{x) = TSpp: . {A3) 



The field operator in the in-state is defined as 



£p 



d+)i 



S+) 



(-)i 



S-h 



ix)af> + (l))r>{x)af>] 



(A4) 



being Hermitian conjugated, i.e., [ap 



{±)i+ _ *W 



The diagonahzation of the Hamiltonian is 



achieved by the transition to quasiparticle representation. The Bogohubov transformation 
for creation and annihilation operators of quasiparticles has the following form: 



^(-) 



b)r>{t) 



^(+) 



b'^Jit) 



with the condition 



= ap(t) ap + (3p{t) a:^p 
= a*p{t)a^^^ + (3;{t)a'p-^ 

\ap{t)\' - \(3p{t)\' = 1 . 



(A5) 



(A6) 



The field operator (A4) may be rewritten now in terms of the quasiparticle operators for 
the creation and annihilation processes 



{X 



d^p[<^f\x% 



<|.W(x)4+)] 



(AT) 



.(±) 



where the new amplitudes $^ (x) are related to the functions (Al) by a canonical 
transformation similar to (A5) 

(f)'f\x) = ap{t) 4-)(x) + (3*p{t) 4+)(x) , 



(A^ 



These transformations should be considered together with Eq.(A5). 

Let us find now the equation of motion for coefficients of the Bogoliubov transforma- 
tion. It is assumed by analogy with Eq.(Al) that 



$f)(x) 



L-3/2Q(±)(p,t)e±^®We^P"^ 



(A9) 



where Q^^\p,t) are new unknown functions and the dynamical phase G(t) is defined by 
Eq.(|l2D. After the substitution of Eq.(A9) into Eq.(A8) we have 



g^~Kp,t) = ap-(t)g(-)(p,t)e-eW + /3|(t) Q(+)(p,t) e^^W , 



(AlO) 



where functions g^^\p^t) should obey Eqs.(A2). According to the Lagrange method (e.g. 
p3[]), this requirement may be associated with the following additional relations 






iuj{t) 
iuj{t) 



c^iit) Q^+\p,t) 



je{t) 



- (3pit) Q^'\p,t) e 



ie{t) 



(All) 



ap (t) Q'--\p,t) e-^""^'^ + /3|(t)Q(+)(p,t)e^®W 



After differentiation of Eqs. (All) and the use of the relations (AlO) and (A2) we obtain 
the equations of motion for the coefficients of the canonical transformation (A5): 

"^^') = H) ^^^'^ '"'''' ' ^'^'^ = ^ "^-^') '"''"' ^^''^ 

and the relation Q^^\p^t) = [2a;(t)]^^/^. Now, introducing new creation and annihilation 
operators by 

cf)(t) = 6f)(t) exp[±ze(t)], (A14) 

and using the Eqs. (A5) and (A12) we get the Heisenb erg-like equation of motion |[T8| , |T9t 

^^^ - -^ -^it) + 1 [H{t), cf\t)] , {AU) 



20. 21 



dt 2uj{t) 

with the Hamiltonian of a quasiparticle system, 

H{t) = J2uj{t) 0^\t) cH(t) + S(_-)(t) c^^J(t)} . (A15) 

p 

The canonical transformation (A5) does not change the form of the commutation relations 

iHr^it), c'^\t)U = [clr\t), c(|)(t)]_ = 6pp,. {A16) 

Combined with Eqs.(Al5) and (A16), the equations of motion (A14) are sufficient for 
deriving the kinetic equation ( p2D with the source term (pG]). The derivation procedure is 
very similar to the case of the electron-positron system considered in Sect. 3. 
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